ON THE CONVERGENCE OF KERGIN AND HAKOPIAN INTERPOLANTS AT 

LEJA SEQUENCES FOR THE DISK 
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ABSTRACT. We prove that Kergin interpolation polynomials and Hakopian interpolation poly- 
nomials at the points of a Leja sequence for the unit disk D of a sufficiently smooth function 
/ in a neighbourhood of D converge uniformly to / on D. Moreover, when / E C°°(D), all the 
derivatives of the interpolation polynomials converge uniformly to the corresponding derivatives 
of/. 

o 

^ ■ 1. Introduction 

C . 

Kergin and Hakopian interpolants were introduced independently about thirty years ago as 
natural multivariate generalizations of univariate Lagrange interpolation. The construction of 
these interpolation polynomials requires the use of points, usually called nodes, with which one 
obtains a number of natural mean value linear forms which provide the interpolation conditions. 
Kergin interpolation polynomials also interpolate in the usual sense, that is, the interpolation 
polynomial and the interpolated function coincide on the set of nodes but this condition no 
longer characterizes it. The general definition is recalled below. Approximation properties 
of Kergin and Hakopian interpolation polynomials have been deeply investigated, see e.g., 
03 |3l SI [5J 19|. Elegant results were in particular obtained in the two-dimensional case when the 
nodes forms a complete set of roots of unity (viewed as a subset of IR 2 ). Thus, in Ifl2ll . Liang 
established a formula for Hakopian interpolation at the roots of unity in M?, and later together 
^ ■ with Lii lfl4l estimated the remainder and proved that Hakopian interpolation polynomials at 
the roots of unity of a function of class C 2 in a neighbourhood of the closed unit disk Del 2 

O 
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converge uniformly to the function on D. Thanks to Liang's formula, further authors investi- 
gated (weighted) mean convergence of Hakopian interpolation (see lQ"4"l[T3l ). On other hand, 
in 1997, using a beautiful formula for Kergin interpolation at nodes in general position in M 2 , 
Bos and Calvi (61 independently established a similar convergence result for Kergin interpola- 
tion. If C n denotes the set of n-th roots of unity, J4?[C n ; ■] (resp. J^[C„; •]) the Hakopian (resp. 
Kergin) projector, the results can be stated as 

JT[C„;/] ->■ / and Jt?[C n ;f] -> /, uniformly on D, for every / G C 2 (D). 



X 



In the above results, going from n to n + 1, we need to change all the nodes and it seems natural 
to look for similar results in which C n would be replaced by a set E n such that E n C E n +\, which 
comes to find sequences of nodes rather than sequences of arrays of nodes. It is the purpose of 
this note to exhibit such sequences. They enable us to obtain series expansions of the form 

/(*) = E V ( e o,---,e rf ,D rf /(-)(x-e ,...,x-e rf _i)) 

for / G C°°(D), where the fJ,(eo, . . . , ej, •) are certain mean value linear forms (whose definition 
will be specified below) and D d f(a) denotes the d-th total derivative of /. The sequences 
that we shall use are Leja sequences for D and the results of the present paper have been 
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made possible by recent progresses on the study of Leja sequences (and associated constants) 
contained in |[2l [H Q. To treat the case of Hakopian interpolation, we shall prove a formula 
for Hakopian interpolation at nodes in general position in R 2 which reduces to Liang's formula 
when the nodes form a complete set of roots of unity and which is of independent interest. The 
proof of our convergence results requires a somewhat higher level of smoothness than in the 
case of interpolation at the roots of unity. The question whether we can weaken the smoothness 
of the interpolated function is still unanswered. Moreover, when the interpolated function is 
in the class C°°, we show that all the derivatives of the interpolation polynomials converge 
uniformly to the corresponding derivatives of the interpolated function. 

Notations. The scalar product of x = (x , .. -,x N ) and y = (y , . . . ,yF) in M. N is defined by 
(x,y) := Y!j=i an d the corresponding norm of x is ||x|| = a/ (x, x) . Let K be a compact set 
in M. N . For each continuous function / on K we set = sup{|/(x)| : x £ K}. The space 
of fc-times continuously differentiable functions on a neighbourhood of K is denoted by C k (K). 
For / eC k (K),k> l,we set 

d\«\f 

D0!/= (dx^ ■ ■ ■ (dx»)^ « = \a\ = a l + --- + a N <k, 



D y f(x) = Df(x)(y) = £ |4 (x)y j , xC K. y ( v 1 v 

Pl dxJ 



The euclidean norm of the linear form D/(x) is denoted by ||D/(x) || . We have 

\\mx)\\ = \t(UiMf 



dxJ 



We also denote by @>d(M. N ) the space of polynomials of variables and degree at most d. 



2. The definition of Kergin and Hakopian interpolants 

It is convenient to recall some definitions and properties of interpolation polynomials in their 
full generality. In particular, we shall introduce Kergin and Hakopian interpolation polynomials 
as particular cases of a more general procedure. 

Given a convex subset SI C M. N and a tuple A of d+ 1 not necessarily distinct points in Q, A = 
(ao,ai, . - £ CL d+1 , the simplex functional / is defined on the space of continuous 

[a ,...,a d ] 

functions C(Q.) by the relation 

J f:= Jf(a +Y,tj(aj-ao))dt, feC(Cl),d>l, (2.1) 

[a ,...,a d ] A d 7 1 

where dt = dti ■ ■ ■ dtj stands for the ordinary Lebesgue measure on the standard simplex = 

{(t\,t2,...,td) e[0, l] rf , Y?i=\tj < !}• In the case (i = we set / / = /(a ). 

M 

The following theorem leads us to the definition of mean-value interpolation. Its proof can 
be found in iflOll or BU. 

Theorem 2.1. Let Q. be an open convex subset of~R. N , A = (ao, . . . ,aj) be a tuple in Q. and 
let k £ {0, . . . , d}. For every function f £ C d ~ k (Ci), the set of all (d — k) -times continuously 
differentiable functions on £1, there exists a unique polynomial P on Mr of degree at most d — k 
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such that 

J D a (f-P) = 0, \a\=j, j = 0,...,d-k. (2.2) 

[a a ,...,aj +k ] 

Definition 2.2. The polynomial P in (12.21) is called the k-th mean-value interpolation polyno- 
mial of / at A and is denoted by if ^ [A; /] or if W [oq, • • • , a d ; f] ■ 

There is an explicit but rather complicated formula for mean-value interpolation, see IfTOl 
Theorem 1] or [[91 Theorem 4.3]. Here we summarize a few basic properties of mean- value 
interpolation. 

(1) The polynomial ifw [A;f] does not depend on the ordering of the points in A, 

(2) The operator if ^ [A] : / e C rf_ *(£2) i — ► if W [A;/] G & d - k (R N ) is a continuous linear 
projector (when C rf_ *(Q) is equipped with its standard topology), 

(3) For each / G C d ' k {Q), the map A G i — ► if W [A;f] G &> d - k (R N ) is continuous, 

(4) For any affine mapping *P : — >• M M and any suitably defined function /, we have 
if W [A; / o = if « [¥(A) ;/]o«P 

(5) The polynomial if (°) [A;/] interpolates / at the a/s and becomes Taylor polynomial of 
/ at a of order J when ao = ■ ■ ■ = a d = a. 

The most interesting mean-value interpolation polynomials are Kergin interpolants which 
correspond to the case k = 0, 

Jtr[ao,...,a d ;f) =Sf®[ao,...,a d ;f\, (2.3) 

and Hakopian interpolants which correspond to the case k = N — 1 and d > N — 1 , 

J^[ao, ...,a d ;f}= Sf<- N ~V [a , ...,a d ;f}. (2.4) 

When the ay's are in general position in M. N - that is, every subset of Af+ 1 points of A defines 
an affine basis of M. N - then Kergin operator extends to functions of class C , see (6l p. 206- 
207]. On the other hand, under the same condition on the points, Hakopian interpolation is 
characterized by the following relation. For P G ^d-N+i (^ N )> 

P = J?[a ,... : a d ;f}^ J (f-P) = 0, < h < ■ ■ ■ < i N < d. (2.5) 

Hence, in that case, derivatives are no longer involved and the Hakopian operator extends to 
continuous functions. 

3. Error formulas for Hakopian and Kergin interpolants in M 2 
We now restrict ourselves to the two-dimensional case. 

3.1. For x = (x ,jt) G M 2 , we denote by x := (— x 2 ,x l ), the image of x under the rotation 
of center the origin and angle n/2. As usual, to x = (x ,jt), we associate the complex number 
x l + ix 2 with i = a/— 1 which we still denote by x. With this notation we have x = ix. Assume 
that the points a,- are in general position (so that no three of them are aligned). We consider an 
one-variable polynomial of degree d — 1 defined by 

d-l 

h st (w) = | | (w- ({a s -a t ) L ,a m )), w G E, s ^ t, < s,t < d- 1. (3.1) 

The polynomial h st appears in the formulas for Kergin and Hakopian interpolation polynomials 
and plays an important role in our arguments. It is worth pointing out that h st is a multiple of 
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the polynomial q ts used in relation (2.18) in [6]| where basic properties of q ts are established. 
To make our exposition self-contained we state and prove a few properties of h st . 

Lemma 3.1. Let d>3 and let A = (ao,ai, . . . ,cid-i) be a d-tuple of points in general position 
in IR 2 . Then 

(1) h st {{(a s -a t ) L ,a v )) = 0, 0< v< d-l; 

(2) ^(((as-at) 1 ,^)) =U d m J 0)m ^ t ((a s -a t ) ± ,(a s -a m )); 

(3) Ifu, v and s are pairwise distinct and ((a s — at) 1 , (a u —a v )) = 0, then 
ti st (({a s -a t ) L ,a u )) =0. 

Proof. Observe that h st (((a s - a t ) L ,a x )) = U^o tm ^ s {(a s ~ at) 1 , (a v - a m )). The product in 
the right hand side has the vanishing factor ((a s — a t ) ,a v — a v ) when v^s and the vanishing 
factor ((a s — at) 1 , a s — a t ) when v = s. Thus h st (((a s — at) 1 , a v )) = 0. For the proof of assertion 
(2), we just compute the derivative of h st , that is 

Kt(w)= £ (w-das-at) 1 ^)). (3.2) 

m=0,m^=s j=0,j^=m,s 

It is easy to see that the vanishing factor ((a s — at) 1 , (a s — a t )) is contained in the product 
]!j=o j^m s(( as ~ a t) 1 ^ ( a s ~ a j)) whenever m^t. Hence, in view of (13.21) . we have 

d-l 

h'Mifls-at) 1 ,^)) = n ((fl ff -flf) ± ,(a s -fly)). (3.3) 

For the last assertion it is enough to verify that flf-n *(( a s ~ a t) 1 , (a u — a;)) = for every 
m 7^ s. To prove this we only notice that the product in the left hand side contains the vanishing 
factor ((a s — a t ) , (a u — a u )) if m — v and the vanishing factor ((a s — at) 1 , [a u — a v )) if m ^ 
v. " □ 

Theorem l3.3l below gives a formula for Hakopian interpolation polynomial in IR 2 . It is similar 
to that of Kergin interpolation polynomial found by Bos and Calvi. Here, we denote by cv(A) 
the convex hull of the set A. 

Theorem 3.2 (Bos and Calvi). Let A = (ao,ai,...,ad-i) be a tuple of d points in general 
position in the plane. Then the (extended) Kergin operator Jt[A] is continuously defined on 
C l (cv(A)) by the formula 



jr[A;f}= d ff(a J )P J + £ P st [ D (at _ flj) x/, 
y=o o<s<t<d-i 



[a s ,ai\ 



where Pj is the real part of the (j + \ )-st fundamental Lagrange polynomial corresponding to 
the complex nodes ao, . . . , a^-i, that is 

Py(*V) = *( fl (xl +* ) ~ am ), (3-4) 

m—Q.my^ j J m 



and 



P (r) _ hstdias-at) 1 ^)) 

\\a s -a t \\ I h st {((a s -a t ) ± ,a s )) 



Theorem 3.3. Let d>2 and let A = (ao, a\ , . . . , i) be a tuple ofd points in general position 
in the plane. Then the ( extended) Hakopian operator Jtf?[A] is continuously defined on C(cv(A)) 
by the formula 

■^[Mf]= £ Qst I f, (3.6) 



0<s<t<d-l 

\a*,at\ 



where 

_ h'Mas-a^.x)) 
Ust[x — 77-777 7- 77- {j- I) 



Proof. Of course, when d = 2 then <2oi = 1 and (13.61) is trivial. Now, suppose that d > 3. Since, 
for all < s < t < d — 1, Q st is a polynomial of degree at most d — 2, the polynomial defined 
in the right hand side of (13.61) belongs to ^d^O^ 2 )- Let us call H this polynomial. Thanks to 



(3.8) 



(3.9) 



(12.51) . we have to show that 

J H= J f, for all < u < v < d- 1. 

[a u ,a v ] [a u ,a v ] 

In view of (13.61) . it suffices to verify that 

J Q s t = S us d vtl for all < s < t < d- 1, < u < v < d- 1, 

[a u ,a v ] 

where 8 is the Kronecker symbol. Looking at (13.71) . we have 
1 

/ & = /e.,,K+wK- fl „))dw 

[a u ,a v ] 

1 

= 77T77 ri / Kt {(( a s - a t ) L ,a u ) + w((a s - at) 1 , (a v - a u )))dw. (3.10) 

h st ({{a s -a t ) ± ,a s }) J 



To deal with the last integral we examine three cases. 

First, we assume that (sj) = (u,v). Then, since ((a s — at) 1 , (a t — a s )) = 0, relation (13.101) gives 

Q st = l. 



[a s ,a t \ 

The second case occurs when (s,t) 7^ (u, v) and ((a s — at) 1 , (a v — a u )) = 0. Then s ^ u and 
s ^ v. Indeed, if, for exemple, s = u, then the relation ((a s — at) 1 , (a v — a s )) = implies that 
a s ,a t and a v are collinear, contrary to the hypothesis. Now, the integral term in (|3.10l) reduces 
to 

Q = Ktidas-at) 1 ^^) _ Q 
^({(as-at^^s)) 

[a u ,a v \ 

where we use Lemma [3TT 3) in the second equality. The last case is when (s,t) 7^ (u,v) and 
((a s — at) 1 , (a v — a u )) 7^ then, calculating the intergral (13.101) . we have 

_ /?, rf (((a y -a f ) J -,a v )) -h st {{{a s - at) 1 ^^) (3 11) 

h' st (((a s -a t ) 1 ,a s ))((a s -a t ) 1 ,(a v -a u )) 

[a u ,a v ] 

Now, Lemma ISTTlT ) follows that the right hand side of (13.1 II) vanishes and the proof is com- 
plete. □ 
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3.2. We use the above formulas to establish multivariate analogues of the classical Lebesgue 
inequality for Lagrange interpolation. For the proof, we refer to fl6] Theorem 1.1] and Ifl4l 
Theorem 5]. 

Lemma 3.4. Let A = (ao,ai, . . be a tuple of d points in general position in the plane, 

and let K ci 2 be a convex compact set containing A. Then for every f G C 1 (K) and every 
Qd-i e «^_i(M 2 ) we have 

d-l 

\\f-Jt[A-,f}\\ K <(l+Y,\\Pj\\K)\\f-Qd-i\\K+diam(K) £ \\P st \\ K \\Df -DQ d ^\\ K , 

j=0 0<s<t<d-l 

(3.12) 

where diam(^) is the diameter ofK, the polynomials Pj and P st are defined in < \3.4\) and di.51) 
respectively. 

Lemma 3.5. Let A = (ao,a\,... , a^_i) be a tuple of d points in general position in the plane, 
and let K cK 2 be a convex compact set containing A. Then for every f G C(K) and every 
Q d -2 e & d -20& 2 ) we have 

\\f-je[A;f]\\ K <(l+ £ ||&»|k)||/-&/-2||jr, (3-13) 

0<s<t<d-l 

where the polynomials Q st are defined in A3. 71) . 

4. Kergin and Hakopian interpolants at Leja sequences for the disk 

Definition 4.1. Let D be the closed unit disk in the complex plane and E = (e n : n G N) be a 
sequence of points in D. One says that £ is a Leja sequence for D if the following property 
hold true, 

d-l d-l 

| Y[ ( e d - ej) | = max | (z — ef) \ , for all d > 1 . 

j=o zeD j=0 

A J-tuple Ed = (eo,e\, . . .,e^-i) is called a J-Leja section. In this paper we only consider 
Leja sequences whose first entry is equal to 1. It is not diffcult to describe the structure of Leja 
sequences for D. The following theorem is proved in [0. 

Theorem 4.2 (Bialas-Ciez and Calvi). The structure of a Leja sequence E = (e n : n e M) for 
the unit disk D with eo = 1 is given by the following rules. 

(1) The underlying set of the 2 n -Leja section E^ consists of the 2 n -th roots of unity 

(2) The 2 n+l -Leja section is (E2»,pE^'), where p is a 2 n -th roots of -1 and Ey) is the 
2 n -Leja section of a Leja sequence E^ = (e„ : n e PJ) for the unit disk with e^ = 1. 

Next, we use Lebesgue-type inequalities for Kergin and Hakopian interpolants along with 
the method of Bos and Calvi to prove the following convergence results. 

Theorem 4.3. Let ^\E d \f\ denote the Kergin interpolation polynomial of f with respect to 
the Leja section E^ = (eo, . . . , ej-i) of a Leja sequence E = (e n : n e N) for D. 

(1) If f G C 4 (D), then Jtf[Ej;f} converges uniformly to f on D as d — > °°; 

(2) Iff G C°°(D), then (jf[E d ;f}) converges uniformly to D^f on D as d — >• °°, for 
every two-dimensional index j5. 

Corollary 4.4. For every f G C™(D), the series 

£ / V d f(-,x-eo,---,x-ed-i) 

% ,...,e d ] 



converges to f uniformly on D. Moreover, the convergence extends to all derivatives. 

Proof. In view of Newton's formula for Kergin interpolation (see [15, Theorem 2]), the (d+ 1)- 
st partial sum of the series is exactly Jf[E d ;f] . □ 

Theorem 4.5. Let Jif[E^;f] denote the Hakopian interpolation polynomial of a function f with 
respect to the Leja section Ed = (eo, ■ ■ ■ , e^-i) of a Leja sequence E = (e n : n G N) for D. 

(1) Iff G C 5 (D), then Jif[Ed',f] converges uniformly to f on D as d — > °°; 

(2) Iff G C°°(D), then (jf?[E d ;f\) converges uniformly to D p f on D as d -»■ °° for 
every two-dimensional index /3. 

Corollary 4.6. For every f G C°°(D), the series 

L £ / D f{-^-e Jlt ... i x-e Jd _ 1 ), 

d=10<ji<j 2 <-<jd-i<d-% J o i 
converges to f uniformly on D. Moreover, the convergence extends to all derivatives. 

Proof. Looking at the formula for Hakopian interpolation (see ifTTl [TOlO . the <i-th partial sum 
of the series is exactly ^[E^f]. □ 

Remark 4.7. We denote by F p the set of functions from {0, . . . , d — 1} to {1,2}. For % G Fj, 
we set 05 (t) = (a,b) with a (resp. the number of times that x takes on the value 1 (resp. the 
value 2) and we write (x — e) T := nf=o \( x ~ e i)%(i)> wn erc (x — e/) T (j) is x l — (resp. x 2 — 3e ; ) 
if t(/) = 1 (resp. x(i) = 2). In particular, (x — e) T is a polynomial of degree J. We have 

D rf /(-,x-e ,...,x-e rf _ 1 )= £D a M/(-)(*-e) T . (4.1) 



The series expansion in Corollary 14.41 can be rewriten as 

/(*)=£ L / D a W/(-)(x- e ) T . (4.2) 

Thus the polynomials (x — e) x can be regarded as a generalization of the classical Newton 
polynomials and the above expansion as a multivariate Newton series expansion. A similar 
observation could be done with Corollary 14.61 

In the rest of this note we always interpolate at Leja sections E d . In order to use the Lebesgue- 
type inequalities given in (13.121) and (13.131) . we need a kind of Jackson theorem that we now 
recall. 

For/GC fc (D), we set 

11/11*:= I K/lk 

|j8|<fe 

6)(f;8):=sup{\f(x)-f(y)\:\\x-y\\<8,x,yeD}, 5)= £ o(D^/;5), 5 > 0. 

The following theorem, proved in [fT6l p. 164], is due to Ragozin. 

Theorem 4.8 (Ragozin). Given f in C k (D) with k > 0, there exists polynomials Qj with 
deg Qd < d such that 

Wf-QdWo < M(k)d- k [d- 1 \\f\\ k + ©(/*>; 1/d)] , 
where M{k) is a positive constant which depends only on k. 
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In O, Bos and Calvi slightly modified a method of Ragozin and used Theorem 14 .81 to prove 
a simultaneously approximate theorem for C 2 functions on D (see [6, Lemma 4.1]). Examining 
their proof, we see that the proof easily extends to C k functions. We state the generalized result 
without proof. 

Lemma 4.9. Given f in C k (D) with k > 1, there exists a sequence of polynomials Qd with 
deg Qd <d such that 



lim^ 1 ||/-2 rf ||fl = and lim d h ~ l \\Df -DQ d \\ D = 0. 

d— d— 

Next, we investigate the growth of Lebesgue-type constants 

d-\ 

£ \\ p j\\D, W p *\\d and £ 

y=o o<s<f<rf-i o<i<f<d-i 

-d-i 



Looking at the formula for Pj in Theorem 13.21 we see that Ey=o ll^jllo * s dominated by the 
Lebesgue constant A(£^) for Lagrange interpolation corresponding d complex nodes ej, < 
j < d — 1, over D. But (81 Corollary 7] tells us that A(£' d ) = 0(Jlog d) as J — >• °°. It gives us 
the estimate 

d-l 

Y \\Pj\\ D = 0(dlogd) asi->«,. (4.3) 

j=o 

The estimates for the remaining Lebesgue-type constants are simple consequences of the fol- 
lowing two theorems that are proved in the last section. Here, in the formulas for h st and h' st , 
we take a , 



ej so that 



d-l 



h st {w)= Yl ( w -{i e s-e t ) L ,e m )), w6l, 

m=0,mj^s 



>4«(< 



.» 



d-l 



ra=0,ra7^i',f 



(4.4) 
(4.5) 



Theorem 4.10. We have \\P st \\d < 2d for alld>2 and < s < t < d — 1, w/?ere 

^(((es-eO ± j JC )) 



jcGD. 



|e 5 -e f | 2 /4 (((cs-cO 1 ,^)) 
Theorem 4.11. We have \\Q s t\\D < 4J 3 for alld>2 and 0<s<t<d—l, where 

^(({e g -e t )\x)) 



Qst(x) 



xeD. 



h> t {((e s -e t ^,e s )y 
Next, we note that D satisfies a Markov inequality, that is, 

<{degp) 2 \\p\\ D: pe 





dp 




dp 




max | 


dx l 


D 


dx 2 


J 



2 ), (4.6) 

see ifTTl . Repeated application of (|4.6I) yields estimates for each partial derivatives, 

WT^pWd < (degp) m \\p\\ D , p G ^(M 2 ). (4.7) 

Proof of Theorem POl Using Lemma I4T91 for / G C k (D) we can find a sequence of polynomials 
<2</-i G ^d_i(IR 2 ), J G N*, and a sequence of positive numbers (£„:«£ N) that converges to 
such that 

max{||/-Q rf _ 1 || D ,||D/-DQ rf _ 1 || D }<^ T , d>\. (4.8) 



But Theorem ETO] gives Zo< s <t<d-i \\Pst\\o < d 2 (d - 1) and (g/J gives 1 +Ljzl \\ p j\\D = 
0(d\ogd). It follows from (CTTIl) that 

\\f-^[E d -J]\\ D <(o{d\o Z d) + 2d\d~\))^<^ (4.9) 

where M is a constant. The right hand side of (14.91 ) tends to as d — > °° when / G C 4 (D), 
i.e., k = 4. This follows the first assertion. To prove the second one with the hypothesis that 
/ G C°°(D), we first observe that 

|1 Jt[E n+1 ;f] - Jt[E n ;f] \\ D < \\f - Jtr[E n+1 ;f] \\ D + ||/ - JtT[E H ;f\ \\ D < M{£n+ k -t ^ ■ 



(4.10) 

Applying Markov's inequality in (14771) for p = Jf[E n+ \\f\ - Jt[E n ;f) G ^„ + i(IR 2 ) we obtain 

|d» (^[E„ +1 ;/] - jr [B „ ; /]) || p < +*X»+ 'J™ , lie*. (4.1D 

Now, we choose k = k(f5) = 2|/3 1 + 6. Then the series 

y M(£ n+l +£ n )(n + l) 2 W 
n=\ n 

converges. This follows the uniform convergence on D of the series 

(JC[E V JY) + £ (jr[E n+V J]-X[E n -f)). 
»=l 

Hence (^[E^f]) converges uniformly on D as d — > °°, for every /3 G N 2 . A classical 

reasoning show that if J^[Ed',f) converges uniformly to / on D and {^[E^f]^ converges 
uniformly on D for every /3 G N 2 , then 

d' 3 (jTf^;/]) D^/, uniformly on D, for every /3 G N 2 . 

□ 

Proof of Theorem FOl Using Theorem 14 .81 for / G C k (D) we can find a sequence of polynomi- 
als Qd-2 ^ d>2, and a sequences of positive numbers (8 n : n G N) that converges 
to such that 

\\f-Qd-2\\ D <f k , d>2. 

From Theorem 14. 1 1 1 we have Y,o<s<t<d-i WQAd < 2d 4 (d-\) < 2d 5 - 1. Hence, in view of 
Lemma [331 we get 

\\f-jr[E d ;f]\\ D <(l+ £ \\Q st \\ D )\\f- Qd _ 2 \\ D <l^_. 

Q<s<t<d-\ a 



This estimate is the same as (14.91) . The conclusions of the theorem now follow by repeating the 
arguments in the proof of Theorem 14.3 1 □ 
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5. Further properties of Leja sequences for the disk 

5.1. Decomposition of Leja sections. Let E = (e n : n G N) be a Leja sequence for D. As 
observed in 0, repeated applications of the rule in Theorem 14.21 show that if d = 2"° + 2' H + 
h 2" r with n > ni > • • • > n r > 0, then 

Ed = (E 2 »o , poEj} 2 „ ) = (E 2 "o , poEyi] , piPo£'^2«o-2"i ) (5.1) 
= • • • = (£ 2 «o , poE^l , pipo^S , • • • , Pr-1 • • • PiPoE^nl ) , (5.2) 

where each consists of a complete set of the 2'^ -roots of unity, arranged in a certain order, 
and pj satisfies p ■ = — 1 for all < j < r — 1 . 

For a sequence of complex numbers Z = (zk '■ k G N) we define Z(j :k) := (zj,Zj+\ : ■ ■ ■ ,Zk), 
of course, Z(0 : it - 1) = Z k . For J = 2"° + 2" 1 H h 2 n ' with n > »i > • ■ • > n r > 0, let us set 

d_ 1 = 0, rf / = 2 n ° + --- + 2 n ;, for0<j'<r. (5.3) 
With this notation, in view of (I5.2I) . we have 

: d - 1) = £ 2 » and E(dj : - 1) = p; • • -po^J?, < j < r- 1. (5.4) 

From now on, we always denote 0„ := arge„ for n > and := argpy for < j < r— 1. Since 
p^ ; = — 1 we may put (pj = ^ lq '^if )n , qj G N. The following lemma is similar to jH Lemma 3]. 

Lemma 5.1. Let E = (e n : n G N) be a Leja sequence for D and d = 2 n ° + 2" 1 H h 2" r wzY/z 

«0 > «i > • • ' > n r > 0- Then 

a) nt:L 1 \z-e m \ = \z 2no -i\; 

(2) n^T 1 |z - e m \ = | (zpo- 1 • • • pjY J+1 - 1 1, < 7 < r - 1; 

(3) it^'U 1^ - = 2 " ;+I ' ^ ^ ^ - L -l<i<r-l. 

Proof. Since ^"i+i^ f orms a complete set of the 2 n J +1 -st roots of unity, the relation in (15.41) gives 



dj+i-l 



n i*- g «i= n iz-py-poci= n i(pr--po) ^-^i = i^p 1 ■••p. / 1 ) 2 "' /+1 



m=dj ^f=F^' +1 ' „aF^ + ^ 

1 eeE 2 n J+l e€E fj+l 



This proves the first two assertions. For the third one, we observe that e\ = pj--- poe' with 



e' eE^. It follows that 



dj+l-i 

j j \e k -e m \= j j | e '- e |=2 n v+i, 



m=dj,m^k ^fj+l 

since the middle term is the modulus of the derivative of z — 1 at e . This completes the 
proof. □ 

5.2. Some trigonometric inequalities. Let T n be the monic Chebyshev polynomial of degree 
n, that is 2"~ 1 r„(cos <p) =cos(n<p). If cos(rc/3) 7^ ±1, then the equation T n (x) = r„(cos/3) has 
n distinct roots: cos(/3 + 2mn/n), m = 0, ...,n — 1. Hence 

n-i 

r„(cos<p) — r„(cos/3) = ]~[[cos<p — cos(/3 +2m7r/n)]. (5.5) 

ra=0 
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Since both sides of (15.51) are continuous functions of /3, relation (15.51) holds true for all j8 £ R. 



Now, the relation in (15.41) implies that, for — 1 < y < r — 1 , 

{0 m : dj <m< d j+ i - 1} = {(po + ■ ■ ■ + (pj + 2nk/2 n J+ l [2k] :0<k< 2 n J+ l - 1}, (5.6) 

where we write a = /3[2tt] if a = /3 mod 2% and the <pfs do not appear when j = —1. Using 
(15.51) for n = 2 n i+ l and /3 = <po H 1- <Pj — V we obtain from (15.61) the following result. 

Lemma 5.2. Let E = (e„ : n e N) Z?e a Le/a sequence for D, 9 n = arge„, and d = 2"° + 2" 1 + 
h 2' 1f vwY/z «o > > " ' > n r > 0- Then for all y/ £ R and — \<j<r—\we have 

\\ [cos 9 - cos ( e„ 7 - yr)] = — ^ — f cos 2 n J+ x (p - cos 2"> +1 (<po + • • • + q>j ~ Y) 

m=dj 2^ 

where the di 's are defined in ( 15.31) . 

Lemma 5.3. //n > 1, < j < n — 1 and sin(/3 + 2 j'ft/n) 7^ 0, dien 

^Icycon-a-/,)! < 2 „- 1|sin(g+2W „ )r mi 

Proof. In view of (15.51) we have 



'■i-r 1 /n „ / m cosfnffl) — cos(n/3) 

[ J |cos<p-cos(/3+2m77;/n)| - ' 1 



2"- 1 [cos<p-cos(/3 +2jn/n)] 



(5.7) 



Set if = j8 +2j%/n, y\ = ( 1/2) (<p + iff) , xp~% = (1/2) (<p — i/a). Then the sum-to-product formula 
for cosines transforms the right hand side of (15.71) into 2~' !+1 | sinni//i sinnl/^l/l simpi sin i//" 2 1 . 
Since sin y = sin yf\ cos i/a 2 — cos yf\ sin 1//2, we obtain after simplication 

j I I cos<p — cos(/3 + 2m7j;/n)| 

m—0,mj^j 

1 / 1 sin « i//"2 1 |sinni//i 



sin cos 1//2 — 


cos y/\ sin 1//2 




sinni//i sinnl//2 


2 m-i 


sin yr\ sin sin y/2 





< 



2"~ 1 1 sin \\f\ V I sin 1//2 1 | sin \\f\ \ 
2n 



2"- 1 |sin(/3+2j^/n)|' 

where we use the classical inequality | sinna| < n\ sina| for a £ R, n £ N in the third line. □ 



Now, for dj < s < dj + i — 1, Lemma 1531 and equation (15.61) imply that 

n [cos(p-cos(0„ 7 + /3)]< „. +1 2 ' 2>lJ+1 , (peR. (5.8) 

m =dj,m^s 2 2J - 1 |sin(0, + /3)| 

In (I5T81) . taking /3 = -MA with 5 < we get the following result. 

Lemma 5.4. Under the same assumptions of Lemma 15.21 z/dj < s < dj+\ — 1 with —1 < j < 
r — 1 and s < ? < d — 1, ?/zen 

j +1 8 +6 2 2 n J +1 

[cos<p-cos(0 m -^— '-)]<— -pr- (PGM. 

m=dj,m^s ^ I sill 2 I 
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5.3. Further results. The following lemma will be used to get lower bounds for the denomi- 
nators of P st and Q st . 

Lemma 5.5. Let E=(e n :ne N) be a Leja sequence for D and d > 2, d = 2 n ° + 2' n H h 2" r 

with no > n\ > ■ ■ ■ > n r > 0. Then 11^=0 m^s \ es ~~ e ' n \ — ^' f or < s < d — 1. 

The proof of Lemma 15.51 requires a purely trigonometric inequality given in the following 
lemma. For the proof we refer the reader to 0. 

Lemma 5.6. Let r > 1 and let n$ > n\ > ■ ■ ■ > n r > be a decreasing sequence of natural 
numbers. If (pj = (2qj + \)7t/2 n j withqj G Z, j = 0,...,r— 1, then 

r-l 

]" 1 |sin2'V+i- 1 (<p-(p <p 7 -)| > (l/2'«)|cos2"°" 1 (p|, <p el. (5.9) 

j=o 

Proof of Lemma 1731 The case r = is trivial, since YI^Lq ]„^ s \ e s — ^m\ = 2 n °. Thus we may 
assume that r > 1 . Notice that, since arg e s = 6 S and arg pj = <pj = (2qj + 1 ) n/2 n i , < j < r — 1 , 

\(e s p Q l ■■ -p k Y k+l -l| = 2|sin2'^'- 1 (0,-( P o (p k )\, 0<k<r-l. (5.10) 

First, suppose that r > 2 and s > 2 ,lQ . Then there exists a unique < j < r — 1 such that 
dj <s < dj + \ — 1, where the rf,'s are defined in (15.31) . We write 

rf-l do-l d j+i~ l r-l d k+1 -l 

Y[ \ e s~e m \= Y[ \ e s~ e m\- Yi \ e s~ e m\- Yi I! \ e s- e m\- (5-1 1) 

m=Q,mj^s m=d-\ m=dj,mj^s k=0,ky^j m=d) c 

We will treat three factors in (15.1 II) independently. The first and the third part of Lemma 15.11 

give rtl, \e s -e m \ = |ef° - 1 1 = 2, and \t^ d ] m+s l e * ~ e ™l = 2 " j+1 ■ 0n tne other hand ' the 
second part of Lemma IBTTl along with equation (15 .101) yields 

r-l r -i 

n n i^-^mi = n \^Po i ---p k i f k+i -M (5-12) 

k=0,kj^ j m=dk k=0,k^ j 

r-l 

= 2 r - 1 |sin2"w- 1 (0,-(po <p k )\. (5.13) 

Since dj < s < d j+ \ - 1, relation (1531) tells us that 5 = <p H h <?>/ + 2qn/2 n j+ l [2n\ for 

some q E Z. To estimate the product in (15.131) we proceed as in (7, Subsection 4.4]. For the 

convenience of reader, we reproduce the proof . For < k < j, we have 6 S — (po <Pk = 

<Pk+i H 1- <Pj + 2qn/2 n i +1 [2n]. Hence the hypotheses on the values of <po, • • • , <p r -i give 

2 n ^-\e s -( Po cp k ) = ((2q k+l + \)7t/2)[7t}. (5.14) 

It follows that 

| sin2'^'- 1 (0, - (po <Pk)\ = 1- (5.15) 

On the other hand, using Lemma [531 for cp = S — (po <P/, i.e., 9 = 2qn/2"j+ l [2n], we 

obtain 

r-l 

YI |sin2' lA+l " 1 (0,-(p o (p k )\ > (\/2 n j +1 - nr )\cos2"j +1 - l (p\ = l/2^+^' v . (5.16) 

k=j+l 
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Combining (15.151) and (15.161) we get 

r-l 

Y[ |sin2"' : + I " 1 (0 s -(po <Pk)\ > l/2 n ' +l ~ nr > l/2 n J +1 . (5.17) 

k=0,k^j 

Note that when j = r— 1, then the left hand side of (15.171) is equal to 1 and inequality (15.171) 
is obviously true. When j = 0, then the factor in (15.151) does not appear. Now, using relation 
(15.171) in (15.131) we get the following estimate 

r-l 

II II \e s -e m \>2 r - l 2- n ^. 

k=0,ky^ j m=dj l 

In this case, we finally obtain 

d-l 

Yl \e s - e m \ = 2 ■ 2 n J+ l ■ 2 r ~ x ■ 2~ n ^ = 2 r . 

We now treat the case r > 1 and < s < 2 n ° — 1 . The proof is the same as above. Indeed, thanks 
to Lemma IBTTl and (15.101) . we can write 

d-l rfo-l r-l 

Y[ \e s -e m \ = j | |ej-emin Y[ \ e ^~ e m\ (5.18) 

m=0,mj^s m=d_\,m^s k=0 m=d] i 

r-l 

= 2"°2'nisin2"*+ 1 - 1 (0,-<p o <p k )\. (5.19) 

k=Q 

Now using Lemma 151)1 in (15.191) we get 

d-l 

]1 \ e « -em\> 2 n »2 r 2- n » +nr \cosT»- l e s \ = 2 r+nr > 2 r , 

m=0,mj^s 

since | cos2"°" 1 5 | = 1 for < s < 2"° - 1. The last case is when r = 1 and 2' ,( » < s < 2"° + 
2"i — 1 . The proof is simple and we omit it. □ 

Remark 5.7. Take d = 2"° + 1 and s = 2 n ° then nJ"=o 1 l e 2"o -e m \ = |ef$ - 1| = | - 1 - 1| = 2. 
Thus the conclusion of Lemma [531 is optimal in some cases. 

6. Proof of Theorems 14 . 1 01 and 14.111 

We continue to use the notation introduced in the previous two sections. Let E = (e n : n € N) 
be a Leja sequence for D and 9 n = arge„, n > 0. We use the formula for the polynomials h st 
and h' st given in (14.41) and (14.51) . We start with the following simple observations 

_ . U — V Bu+Oy - , . U — G v e u +e V j 

e u — e v = 2ism — - — e~^~ and \e u — e v ) =—2 sin — - — e~^~ for uf^v. 

Set a uv := e~^~ 1 and p uv := —2 sin " 2 '' . Then we immediately see that 

{e u -e v )^ = a uv f} uv and \e u -e v \ = \(e u -e v ) L \ = |j8 MV | for w ^ v. (6.1) 

Lemma 6.1. Lef £ = (e n : n e N) be a Leja sequence for D and d > 2, d = 2"° + 2" 1 H h 2"'' 

wzY/z no > fti > • • • > n r > 0. Then 

2 2r \R | d - 4 

|^((( gj - ef )-L, gj ))|>^^_, 0<5< t <J-l. (6.2) 
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Proof. Thanks to (16.11) we may write 



d-l 



K(((e s -e t )\e s ))\ = [J \({e s - e t ) L ,{e s - e m ))\ ( see 63 

m=0,my^s,t 

d-l I „ \ d-l 



= IP* 



n 

m=Q,m^s,t 



a 



i"- 2 n (««.^> n 

m=0,m^s,t \ " m / m =0,m^s,t 

d-l I „ \ d-l 

\d-3 



'St") ' 



-ft 



n 

m=0,m^s 



(6.3) 



e _ e K6s+6m) i{x+e s +e m ) Ws+Ot) 

Since S _ R m — le 3 = e 2 and a st = e 2 t we nave 



a 



Stl ' 



cos( 



71 + S + 6 m S + t , 



sin- 



Q t — G n 



\e t -e„ 



(6.4) 



Combining (16.31) and (16.41) we obtain 



d-l 



\h' st (({e s -e t )^,e s ))\ = \^ \[ 

m=0,m^s,t 
\P Sl \ d - 4 ^ 



\e t -e ri 



n w 



> 



2d-2 

2 2r \M d ~ 4 



- n 

m=0,m^s 
d-l 

n 



2 d 



-2 



where we use Lemma 1531 in the third line. 

Lemma 6.2. Under the same assumptions of Lemma 1(5.7 1 we have 



\h st (({e s -e t )^,x))\< 



2 2r+l d\j3 st 

2 d-2 



d-2 



-, < s< t < d-l, \\x\\ < 1. 



□ 



(6.5) 



Proof. Let us set (a st ,x) = cos (p. Since (e s — e t ) L = ttstfist and (tt st ,e m ) = cos(0 m — G '\ Qf ), in 
view of (14.41) . we get 



\h st (({e s -e t )^,x))\ = |A 



1 d-l 



d-l 



|(a rt ,(x-e m ))| 

d-l 

\Pst\ d ~ l 11 |cos<p-cos(0 m 



e s +e t 



(6.6) 



There exists a unique —1 < j < r— 1 such that dj < s < dj + \ — 1, where the J,'s are defined in 
(15.31) . We can write the trigonometric expression in the right hand side of (16.61) as follows 



r-l 4+1-1 

n n 

k=-l,k^j m=d k 



cos<p — cos(0„ 



n"' 

m=dj,m^s 



cos<p — cos(0„ 



(6.7) 



Thanks to Lemmas 15.21 and 15.41 the first factor and the second factor in (16.71) are dominated 
respectively by 

r ~ x 2 2-2 n J+ l 



n -. 2"* 1-1 

k=-i,tyj z 



and 
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22"^-i| sin ^| 



Combining these estimates with (16.61) we obtain 

r-\ 9 9«/+i ? 2r + l rl\R 



-2 



112^-1 |sin^|- 2<* 
here we use the facts that |j8 sf | = 2| sin | and 2 n ;'+ I < d. □ 
Lemma 6.3. Under the same assumptions of Lemma 1(5. i I we /zave 

|^(((^- ef ) ± ,x})|< ^T 1 , 0<s<t<d-l,\\x\\<l. (6.8) 

Proof. Let £ : M 2 — >■ IR be the linear form defined by £(x) := ((e s — e t ) ,x). Since \(e s — e t ) \ = 
\p a \, we have £(D) = [-|j8 rf |, |jS rf |]. It follows that 

sup|/i, f (((e,-e f )- L ,x))| = || [-1^1,1^0 and sup |/v (((e, - e*)^*)) I = II^4IIhW,IM- 
Since deg/z. rt = d — 1, classical Markov's inequality gives 

KHh/UIW] - ( l /\Pst\)(d- l) 2 ll^f ll[-|ft,|,|fe|]- 
Hence Lemma l6\2l yields 

(H—\\ 2 ? 2r+1 J 3 lft \ d - 3 

sup|/4«(^-e f )\*))l < ^—±L S u P \h st (((e s -e t ^,x))\ < ' ! 



|/3,| -£>'-VU- ^ 2 ,_ 2 



□ 

Proof of Theorems \4. 1 0\ and \4.11\ Lemmas I6TT1 and I6T21 give an upper bound - — 4^ — for the 



2</-2 

numerator of P st {x) and a lower bound — — for the denominator of P st {x). Thus \P st (x)\ < 
2d for all x E D and 0<s<t<d — 1. At the same time, Lemmas 16.11 and 16.31 follow that 
\Q st (x)\ < 2d 3 \j3 st \ < 4-d 3 for allxeD and 0<s<t<d — \. □ 
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